Scaling law for seismic hazard after a main shock 
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After a large earthquake, the hkehhood of successive strong aftershocks needs to be estimated. 
Exploiting similarities with critical phenomena, we introduce a scaling law for the decay in time 
following a main shock of the expected number of aftershocks greater than a certain magnitude. 
Empirical results that support our scaling hypothesis are obtained from analyzing the record of 
earthquakes in California. The proposed form unifies the well-known Omori and Gutenberg-Richter 
laws of seismicity, together with other phenomenological observations. Our results substantially 
modify presently employed estimates and may lead to an improved assessment of seismic hazard 
after a large earthquake. 
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Earthquakes tend to occur in clusters A large 

earthquake often precedes many further aftershocks, 
some of which may themselves be consequential. For in- 
stance, the Northridge earthquake in Southern California 
set off markedly increased seismic activity in the area, as 
apparent in the time series shown in Figure Q In popu- 
lated regions like Southern California that are prone to 
large magnitude events, the likelihood of occurrence of 
strong aftershocks requires accurate assessment in order 
to determine the risk of further damages 0, 0, ■ 

Earthquakes arise from a complicated nonlinear me- 
chanics Q M and individual events are not pre- 
dictable P, Yet, siinple scaling laws describe their 
statistical properties P, S • Empirically, the number 
n{t) of aftershocks occurring at time t after a large maj 
nitude event obeys the modified Omori law 
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The constants K and c are positive, and the exponent p is 
usually found to be close to one. The Gutenberg-Richter 
(GR) distribution 14] for the frequency of earthquakes 
with magnitude greater than m in a seismic region is 



Nirn) = A 10 



— bm 



(2) 



where A is a constant and the critical exponent b typ- 
ically takes a value close to one. Since the magnitude 
is proportional to the logarithm of the size of the earth- 
quake, Eq. (|2Jl describes a scale-free distribution. 

Comparisons between seismicity and critical phenom- 
ena in other physical systems 0, 0, 0| suggest that 
these laws, usually considered to be independent, should 
connect inextricably together. Progress toward an orga- 
nizing theory of seismicity could therefore emerge along 
lines similar to those followed in the modern approach to 
critical phenomena 0|. Here, we show that, in fact, the 



rate of earthquakes following a main shock obeys a uni- 
fied scaling law. This law relates the expected number 
of aftershocks greater than a certain magnitude to the 
magnitude of the main shock itself, and to the time since 
the main shock. It interpolates between the Omori law 
at short times and the GR relation at late times, pro- 
viding a coherent framework from which several other 
phenomenological observations can be deduced. In addi- 
tion, the mathematical form for the rate of aftershocks 
differs significantly from the presently employed formula 
used to assess seismic hazard y, Q . 

Analysis of earthquake recurrence times use a mag- 
nitude threshold above which events are registered |l,3|. 
This magnitude threshold always exists, even if not ex- 
plicitly specified, due to limitations in reliably detecting 
all small earthquakes in a seismic region. Recently, Bak 
et al considered the magnitude threshold itself as a vari- 
able entering into a scaling law for waiting times between 
subsequent events 10]. As we show below, threshold vari- 
ables also enter into a scaling theory for rates of after- 
shocks. 

In order to quantify increased seismicity following a 
large earthquake, we separate smaller events from larger 
ones by imposing two different magnitude thresholds, as 
indicated in Figure The threshold M defines the large 
events. Earthquakes larger than magnitude M are re- 
ferred to as main shocks. A smaller threshold, m, selects 
the remaining events in the catalog to be counted. These 
remaining events, larger than magnitude m, are referred 
to as aftershocks. The quantity RM.m{t)dt is the average 
number of events in the catalog with magnitude larger 
than m occurring in the interval {t, t + dt), given that an 
isolated main shock of magnitude larger than M occurred 
at t = 0. 

We analyzed the earthquake catalog of Southern Cal- 
ifornia from 1984 to 2002 in order to measure the 
rate i?M,m(i)- This catalog includes more than 3.5 x 10^ 
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earthquakes. We identified all events larger than magni- 
tude M (5.5 < M < 6.5). For each main shock, we mea- 
sured, as a function of time, the rate of subsequent events 
larger than magnitude m with 2.5 < m < M — 1.5. The 
resulting data therefore contain background seismicity, 
unlike some other investigations 0|. The measurement 
of the rate following a given main shock was stopped 
when another event of magnitude M or larger occurred. 
This ensures that each aftershock sequence does not in- 
clude events from other main shocks. Some aftershock se- 
quences, therefore, comprise more events and last longer 
than others. To avoid the effects of correlations between 
main shocks, we considered only isolated ones. We ex- 
cluded all events that follow the second of two successive 
main shocks, where both main shocks occur within a time 
interval (of the order of 4 months) of each other. An av- 
erage over main shocks was then computed for each pair 
(M, m) to obtain RM,m{t). 

Depending on the value of M, the same earthquake 
becomes a main shock in one rate measurement, and an 
aftershock in another. For example, a magnitude 5.8 
earthquake could be an aftershock of a magnitude 6.7 
earthquake if M = 6, but is treated as a main shock 
if M = 5.5. This approach reflects the view that af- 
tershocks do not differ from other earthquakes and can 
therefore lead to aftershocks themselves 20] . In this way, 
the pattern of seismicity in space and time can be con- 
sidered to be a single, hierarchically organized, critical 
process [lolll5|. 

Figure El shows the empirical results. At short times 
the rate, i?M,m(i), approaches a universal constant value, 
Rq ~ 10~^, consistent with the presence of the parameter 
c in Omori's Law, Eq. (Q. Strikingly, the rate at short 
times, _Ro: appears to be independent of both M and m, 
indicating that soon after a main shock earthquakes of all 
magnitudes are equally probable. This observation con- 
trasts with the rate estimate used in Refs. , where 
the rate at short times depends on both thresholds M 
and m. However, the constant regime persists longer for 
increasing M and/or for decreasing m. Later in time, a 
1 /t power law decay sets in. Each rate eventually reaches 
a stationary value at late times. The stationary rates for 
different thresholds m agree with the GR law, and do not 
depend on M . The duration of aftershock activity, t^, is 
defined as the time for the aftershock rate, RM,m{t), to 
decline to its stationary level. This level is referred to as 
background seismicity. Comparing Figs.Et^ andEb shows 
that the duration of aftershock activity increases as the 
magnitude threshold of the main shock M increases. 

All these observations can be combined into a single 
scaling hypothesis for RM,m{t) 

^m.»(^)=( ^^^^q,'(„,_,,)^ )g(A10-^^^0 . (3) 

For RM,m{t) to approach a stationary value at late times, 
the scaling function G{x) must behave as G{x) cx x for 



large x. In that case the stationary rate at late time is 
proportional to 10"'"", consistent with the GR law. At 
early times, G{x) for small x is evaluated. If this is con- 
stant, the first term on the right hand side of Eq. Q con- 
trols the behavior. This term corresponds to an Omori 
law with p = I. The parameter A and the critical expo- 
nent b were determined from the GR relation, Eq. (0), to 
he 10-2±0-2sec-i and b = 0.95 ± 0.1 for the South- 
ern California data set. The only remaining parameter 
in Eq. Q is the time to. 

The data collapse method can be used to test and po- 
tentially falsify Eq. |31 This is accomplished by rescal- 
ing the time, t by (>1 lO^''^^) and the rate, RM,m{t) by 
(to + 10''(™~^'^) t), to get a dimcnsionless scaling plot for 
G{x). Using to — 6 sec, we find that all the curves col- 
lapse onto a single curve, as shown in Figure El within 
statistical error. The data collapse verifies our scaling 
hypothesis. The value of to has an uncertainty; we es- 
timate 4 sec < to < 12 sec from changes in the quality 
of the data collapse using different values for to. Evi- 
dently, G{x) changes behavior at a turning point x — xo. 
Two regimes are distinguished: a transient Omori limit 
at short times, G(x xo) —>■ const., and a stationary 
GR limit at late times, G{x ^ xo) x. 

According to our scaling hypothesis, the duration of 
correlated aftershock activity depends on the threshold 
magnitude, M, of the main shock as 



As Xo — 10~^ (see Figure EJ this implies, for example, 
that, on average, the excess rate of aftershocks persists 
for about two months following an earthquake of magni- 
tude larger than six. 

The parameter c in Omori's law is a controversial quan- 
tity in the earthquake literature |l3| . Reported values 
range from less than 0.01 days to over 1 day. Previous 
investigations usually analyzed aftershock sequences fol- 
lowing specific main shocks of interest. It was found that 
the value of c decreases with increasing threshold magni- 
tude m of aftershocks included in the analysis. Compar- 
ing Eq. 1^ with the scaling hypothesis lO, we obtain 

c = to 10''(*^-'") . (5) 

This sheds light on the variability of the c- value found in 
previous investigations. 

The value of the time, to, can be related to the di- 
mcnsionless variable, Xo, and the universal short time 
rate, Rq. Indeed the 1/t decay of the rate at inter- 
mediate times implies that i?o c — RM,m{td)td- As 
RM,m(td) ^ AID"'"" and from expressions Q and lO), 
it follows that 

to — ^ — 10 sec , (6) 
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which is within the range of uncertainty previously esti- 
mated. This provides a consistency check on our statis- 
tical analysis. Conversely, the scaling behavior, Eq. jSJ), 
of the c-value arises from the condition that the rate at 
short times Rq is approximately independent of M and 
m. 

After a large earthquake, forecasts of seismic hazard 
provided by government agencies are based on estimates 
of the rate A m,™ (t) of aftershocks with magnitude greater 
than TO occurring at time t after a main shock of mag- 
nitude M 0,0,0 ■ Within our ansatz, the rate XM,m{t) 
is 

where Rm{oo) is the background seismicity, generally not 
included in aftershock probability evaluation. Eq. jT)) is 
only valid when to < M. It is important to underline 
that \M,m,{t) as derived here in Eq. {Tj), from an em- 
pirically tested law, Eq. (3), differs in several respects 
from presently employed estimates |M |^ Q . For exam- 
ple, it predicts that the rate of aftershocks at early times 
does not depend on the magnitude of the main shock or 
on the aftershock threshold (as discussed above, in this 
limit '^^gj^ — 0). On the contrary, in the stochastic 
model employed in HQ, XM^t) ~ lO''^*^-") immedi- 
ately after the mainshock. 

In conclusion, wc have introduced and tested against 
the record of recent earthquakes in California, a new 
scaling law that unifies the Omori law for aftershocks 
and the GR relation, previously considered to be inde- 
pendent. Our finding indicates a theoretical framework 
(scaling theory) within which the problem of earthquake 
occurrence should be considered. It provides an empirical 
characterization of earthquake statistics that will be use- 
ful to test theories and models in the future. Our results 
may also have practical implications, namely improved 
hazard assessment following a large earthquake. 
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FIG. 1: Time series of earthquakes in Southern California 
in the month of January 1994. The Northridge earthquake 
of magnitude 6.7 occurred on January 17. The dashed line 
represents a main shock threshold Ad = 6 and the dotted lines 
represent aftershock thresholds m = 3, 3.5 and 4. 
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FIG. 2: Average rate of aftershocks at time t after a main 
shock in Southern Cahfornia. In panel (a) the main shock 
has magnitude greater than M = 5.5, in panel (b) greater 
than M = 6.5. Averages are made over (a) 29 and (b) 4 main 
shocks. Symbols corresponds to different threshold magni- 
tudes for aftershocks, respectively (•) m = 2.5, (■) m = 3, 
(♦) m = 3.5, (A) m = 4, (+) m = 4.5 and (*) m = 5. The 
exponent of the straight line is —1. 
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FIG. 3: Data collapse for the rate of aftershocks following a 
main shock in Southern California. This defines the crossover 
function G{x). The numerical values of the parameters are 
b = 0.95, A — 10~^ sec~^ and to = 6 sec. Note that only to is 
a fitting parameter, since the others are determined from the 
GR relation. Symbols are as in fig. 2, colored in black, red 
and green respectively for main shocks of magnitude greater 
than M = 5.5, M = 6.0 and M = 6.5. The exponent of the 
straight line is 1. 



